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ABSTRACT 

The average age difference between pairs of star clusters in the Large MageUanic 
Clouds increases with their separation as the ~ 0.35 power. This suggests that star 
formation is hierarchical in space and time. Small regions form stars quickly and large 
regions, which often contain the small regions, form stars over a longer period. A 
similar result found previously for Cepheid variables is statistically less certain than 
the cluster result. 

Key words: open clusters and associations: general — Magellanic Clouds — stars: 
formation — turbulence 



1 INTRODUCTION 

A previous study of the positions and ages of Cepheid vari- 
ables in the Large Magellanic Cloud (Elmegreen & Efre- 
mov 1996) found that Cepheids closer to each other were 
also more likely to have the same age, as if star formation 
proceeds in small regions faster than in large regions. Be- 
cause Cepheids are relatively old (> 27 Myr in the LMC 
sample, with most around 100 Myr) and could have drifted 
slightly from their point of origin, and because Cepheid ages 
are rather uncertain, we repeat this LMC study here using 
star clusters. The Bica et al. (1996) catalog of clusters in 
the LMC contains 590 members with accurate positions and 
UBV colors, and these colors have been converted to age by 
Girardi et al. (1995). Thus we can use the Bica et al. clusters 
to study the hierarchical properties of star formation in the 
LMC. 



2 CLUSTER AGE DIFFERENCE VERSUS 
SEPARATION 

We consider all pairs of clusters in the LMC for clusters 
within certain age ranges, and determine the average age 
difference At among these clusters as a function of their 
separation 5* for regular intervals of separation. Separation 
is defined to be the deprojected distance between the two 
clusters, considering zero depth to the LMC and an inclina- 
tion of 33° (Luk & Rohlfs 1992). For the age difference, the 
absolute value is used to avoid negative numbers. 

The results are shown in Figure ^. Four age intervals 
are considered: 1 to 100 Myr, 1 to 10 Myr, 10 to 100 Myr, 
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Figure 1. The average age difi^erences between pairs of star clus- 
ters are plotted versus their deprojected angular separations for 
clusters in the LMC. Four intervals of cluster ages are considered. 
For the larger age intervals, the clusters show a correlation in 
the sense that close clusters have more similar ages than distant 
clusters. 



and 1 to 1000 Myr. The average age difference between pairs 
of clusters increases systematically with their spatial sepa- 



© 1998 RAS 



2 Y. Efremov and B. Elmegreen 



1 My < age < 100 My 
Scrambled ages 



10 My < age < 100 My 
Scrambled ages 



1 My < age 
Scrambled b 



10 My 




1 My < age < 1 000 My 
Scrambled ages 



0.01 0.02 0,04 0,1 0.2 0.4 

Separation in Degrees 



0.01 0.02 0,04 0,1 0.2 0.4 

Separation in Degrees 



Figure 2. The same as in figure 1 but with cluster ages scrambled 
randomly among all of the clusters in the sample. The correlation 
disappears for randomized ages, suggesting that the trend found 
for the real data is statistically significant. 



ration. The number of clusters within these age intervals is 
337, 93, 244, and 526, respectively. 

Figure ^ shows a similar study using the same 590 clus- 
ter positions, but now with ages that are mixed up ran- 
domly in order to check purely statistical effects (these are 
not purely random ages, but the real cluster ages that are 
reassigned randomly to different clusters). The correlations 
disappear for this random sample, i.e., the average age dif- 
ference is the same regardless of separation. 

Figure ^ gives the results of 100 random trials using the 
real cluster positions and randomly mixed ages for ages in 
the ranges 1-100 Myr (top) and 10-100 Myr (bottom). Each 
point on the left gives the least-squares fit to the slope of the 
At — S correlation (evaluated in the separation interval from 
0.01 to 1 degree), and each point on the right gives the cor- 
relation coefficient measuring the goodness of fit (l=perfect 
fit). The dashed lines represent the slopes and correlation 
coefficients for the real data, where the clusters have their 
real ages, as shown in Figure 1. 

Figure ^ indicates that the slopes of the At — S correla- 
tions for the real data in the 1-100 Myr and 10-100 Myr age 
intervals differ from random noise in a statistically signifi- 
cant way. The slope is much higher than the slopes in the 
random samples, and the correlation coefficient for the real 
data is good. 

The least-squares fits in the 5* = 0.01 — 1 degree sep- 
aration range for the significant At — S relations shown in 
figure ^ are: 



logAt^yr) = 7.48 -f 0.33 log ^(deg.), 
logAt(yr) = 7.49 + 0.38 log ^(deg.), 
logAt(yr) = 8.20 + 0.42 log ^(deg.). 



(1 - 100 Myr) (1) 
(10 -100 Myr) (2) 
(1 - 1000 Myr) (3) 



The correlation coefficients for these three fits are 0.88, 0.83, 
and 0.82, respectively; the numbers of cluster pairs are 5509, 
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Figure 3. Left: The slope for the least squares fit to the regres- 
sion of average age difference versus cluster separation is shown 
as a point for each of 100 trials with randomly scrambled age data 
and real cluster positions. The slope has an average value of zero 
because there is no At — S correlation for the randomized data. 
The slope for the real data, indicated by a dashed line, is clearly 
larger than the rms slopes for the random data. Right: The corre- 
lation coefficients for the regressions made from randomized data 
are shown as points, and the correlation coefficient for the real 
data is a dashed line. The correlation coefficient for the real data 
is large, indicating that the cluster At — S correlation is statisti- 
cally significant. Two intervals of cluster age are considered. 



2849, and 10614. The 1-10 Myr data has too few clusters 
(470 pairs) to be statistically significant (correlation coeffi- 
cient: -0.0046). 

We conclude from these figures that clusters in the LMC 
form in a hierarchical sequence in which the duration of star 
formation in a region scales with the ~ 0.35 power of the 
region size over scales ranging from at least 15 pc (0.02°) to 
780 pc (1°) at a distance of 45 kpc (Berdnikov, Vozyakova, 
& Dambis 1996; Efremov 1997; Efremov, Schilbach, & Zin- 
necker 1997; Fernley et al. 1998). For the 1-100 Myr age 
range in this 0.01 — 1 degree separation range, this correla- 
tion is 



At(Myr) ~ 3.3S'(pc) 



(4) 



3 AGE DIFFERENCE VERSUS SEPARATION 
FOR CEPHEID VARIABLES 

A similar analysis is now applied to Cepheids in the LMC, 
as in Elmegreen & Efremov (1996), but now we use an age 
calibration that is consistent with the ages of the clusters. 
The ages of Cepheids come from their periods, because both 
the period and the age of a Cepheid correlate with stellar 
mass. The original period-age relation came from Efremov 
(1978) and did not consider overshooting in the stellar evo- 
lution models. The ages of the clusters used in the previous 
section do consider overshooting, so we have to recalibrate 
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Figure 4. Cepheid periods P versus the ages t of tlie clusters in 
which the Cepheids are found, for Cepheids in the LMC. 



the period-age relation using the Girardi et al. (1995) cluster 
age data. 

Table 1 gives the designations, UBV data, and ages of 
select clusters in the LMC, and the names and periods of 
Cepheid variables that appear to be associated with these 
clusters. The UBV data are from Bica et al. (1996), and the 
ages come from the S values as determined in Girardi et 
al. (1995). The sources of the data on the cluster-Cepheid 
associations are indicated by the references. These are the 
most certain associations in the LMC. Many other Cepheids 
are reasonably close to clusters too (Efremov 1978, 1989), 
but their memberships or ages were not judged to be as 
certain. Sometimes it was difficult to select unambiguous 
associations in the rich fields, as was the case for the NGC 
1850 binary cluster and a half dozen Cepheids close to it. 

We have not included the Cepheids with small ampli- 
tudes and sinusoidal curves (subtype Cs), which are either 
the first-overtone pulsators or those at the first crossing of 
the instability strip (Efremov 1968). Their periods do not 
correspond to ages in the same way as for the common 
Cepheids. 

Most of the Cepheids in the table are well inside cluster 
cores, whereas those with references to Efremov (1978) are 
within 2-3 radii of cores. The latter is presumably the case 
for the Alcock et al. (1995) data as well. 

The data in table 1 gives the new period P - age t 
relation shown in figure W. The correlation is 



log t = 8.492 - 0.509 log P, 



(5) 



for t in years and P in days. The correlation coefficient is 
0.90. The result for paper 1 was logt = 8.157 - 0.677 log P, 
which gives shorter ages by a factor of ~ 3. 

The age differences between all pairs of Cepheids in the 
Artyukhina et al. (1995) catalog with ages less than 100 
Myr were determined for intervals of deprojected separation 
equal to 0.01 degree, and the average age differences for 
each interval were found. Figure |^ shows the average age 
difference versus separation for these Cepheids. The figure 
suggests there is no significant correlation between At and S 
in the Cepheid data, but the statistical uncertainty is large 
(out of 1200 Cepheids in the Artyukhina et al. catalog, there 
are only 167 younger than 100 Myr, and among these, there 
are only 1357 pairs with separations less than 1°). 



Figure 5. Average age differences versus separation for Cepheids 
in the LMC. There is no obvious correlation as there is for clusters, 
presumably because Cepheid stars drift from their points of origin 
more rapidly than clusters, and because the Cepheid ages are less 
accurate. 




10 20 30 40 50 60 70 80 90 100 

Trial Number 



10 20 30 40 50 60 70 80 90 100 

Trial Number 



Figure 6. The slopes (left) and correlation coefficients (right) for 
randomly scrambled age data in the Cepheid sample are shown 
as points, while the measured slope and correlation coefBcient are 
shown as dashed lines. The correlation coefficient for the real data 
is as poor as it is for a random sample, and the slope of the At — S 
regression for the real data is within the range of the noise for the 
random data. 



Figure ^ shows 100 random trials for the Cepheids, us- 
ing the observed Cepheid positions with randomly scram- 
bled ages, as for the clusters. The slopes of the At ~ S re- 
lations are shown on the left and the correlation coefficients 
are shown on the right. The dashed lines are the measured 
slope and correlation coefficient, using the real Cepheid data. 
Evidently, the observed relation for Cepheids is consistent 
with noise. If there is a At — S relation for these stars at 
birth, as appears to be the case for clusters, then this re- 
lation is apparently destroyed or obscured by the time the 
stars become Cepheids. 

Possible reasons for the elimination of an initial At — S 
relation include random stellar motions and inaccuracies 
in the Cepheid period-age relation. The correlations shown 
here are better for clusters, perhaps because clusters have 
smaller space velocities than individual stars, so they drift 
less from their points of origin than Cepheids. Smaller space 
velocities for clusters might be reasonable if clusters, born 
in cloud cores, start with the centroid velocities of the star- 
forming clouds, while individual stars, dispersed in associa- 
tions, start with these centroid velocities plus the dispersal 
speeds. 
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4 DISCUSSION 

The positions and ages of clusters in the LMC suggest that 
there is a correlation between the duration of star formation 
in a region and the size of the region. This correlation may 
be the result of star formation inside cloud complexes that 
always live for several turbulent crossing times. 

The turbulent crossing time in a region of star formation 
may be estimated from the size-linewidth relation for the 
clouds that form stars. A compilation of the size-linewidth 
relations for molecular clouds in the Milky Way is shown in 
figure ^ The size S is the FWHM of clouds and clumps in 
various surveys, and the linewidth c is the Gaussian disper- 
sion. The average relation for all the surveys is 

c{kms-^) ^0.7 Sipc)°-\ (6) 

The ratio of 5 to c is shown on the bottom of the figure. 
Half of this ratio gives the turbulent crossing time, 

t,ross^^s{Myr) ~ "j^^^^,^), ~ 0.7S(pc)° '. (7) 
cykm a ) 

Because of the size-linewidth relation, the turbulent crossing 
time increases with cloud size approximately as the square 
root. This slope is comparable to, although slightly larger 
than, the At — S relation for clusters, suggesting that the 
LMC cluster relation results in part from turbulence. For 
separations of 10 pc and 100 pc, the average age difference 
between clusters in the LMC is 3.2, and 2.1 times the cross- 
ing time of galactic molecular clouds on the same scale. This 
suggests that star formation always proceeds with a time 
scale of ~ 2.5 crossing times in the progenitor cloud. 

The size-linewidth relation for clouds in the LMC is not 
well known, but a recent study suggests that the linewidths 
in the L48 clumps are slightly larger for the clump size than 
the hnewidths in galactic GMCs (Kutner et al. 1998). This 
could be a region of unusually high pressure and not normal 
for star formation in the LMC, but if this result is repre- 
sentative, then tcrossingiS) would be slightly lower than for 
Galactic clouds, and the ratio of the star formation time to 
the crossing time slightly larger than ~ 2.5. 

The At — S relation for clusters is not the result of OB 
association expansion or stellar drift at a fixed initial ve- 
locity because then the slope of the correlation would be 1 
instead of 0.3. It has to result from stellar or gas kine- 
matics with velocity dispersions that scale with a fractional 
power of the size of the region. 

There are several important implications for star for- 
mation of the cluster At — S relation: 

1. On all scales over which the correlation exists, small 
regions come and go while large regions continue to form 
stars. This implies there might be some recycling of small- 
scale star-forming material during the lifetime of the larger 
region. Then metal enrichment from supernovae can be 
greater in the most active clouds, and more variable from 
cloud to cloud than previously thought (Elmegreen 1998), 
and total cloud ages can be longer than previously deter- 
mined from cloud disruption times following only one epoch 
of star formation (Leisawitz et al. 1989). 

2. Larger star-forming regions have both larger veloc- 
ity dispersions and larger average ages than smaller regions. 
This trend is similar to that found for OB subgroups and 
whole OB associations and may contribute to the impression 
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Figure 7. Size and linewidth data from molecular cloud surveys. 
Symbols are for giant molecular cloud surveys = dots: Solomon 
et al. (1987), open circles: Dame et al. (1986); quiescent clouds = 
filled triangles: Falgarone et al. (1992), open triangles: Williams et 
al. (1994; the Maddalena-Thaddeus cloud), inverted open trian- 
gles: Lemme et al. (1995; L1498), inverted filled triangles: Loren 
(1989; Ophiuchus), OB associations = crosses: Williams et al. 
(1994; Rosette), plus signs: Stutzki & Glisten (1990; M17). 



that OB subgroups expand into OB associations. There may 
not be this much expansion, however. Instead, there could 
be a difference in the sizes and velocity dispersions of the two 
types of regions from birth. OB subgroups are born small 
and they may stay moderately small during the formation 
time of the other subgroups. All of the subgroups together 
define the association, which is a composite of clumpy sub- 
parts. 

3. The largest regions of star formation in giant spiral 
galaxies, regions like Gould's Belt and other Star Complexes 
measuring 300 pc to 1 kpc, take so long to form (30 Myr) 
that they are not particularly bright on average. They are 
also not unified in appearance by a single bright HIT re- 
gion or concentration of O stars, because the general popu- 
lation of O stars that formed there has already aged off the 
main sequence. Instead, they are visible primarily as con- 
centrations of Cepheid variables and other supergiant stars, 
which is how they were originally discovered (Efremov 1979, 
1989). Most of the O stars are visible only in smaller con- 
centrations, which appear as multiple cores inside the star 
complexes. Thus the largest regions of star formation in gi- 
ant spiral galaxies are often overlooked, especially in Ha or 
UV studies. 

This situation changes in smaller galaxies, where the 
star formation length and time scales are generally shorter 
than in large galaxies (Elmegreen et al. 1996). In small 
galaxies, the largest regions of star formation can form so 
quickly that there are still many OB stars, and then they 
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Figure 8. Schematic diagram showing the size dependence of the 
duration of star formation in various regions. 



appear very bright, like 30 Dor (see review in Elmegreen & 
Efremov 1998). 

4. Regions of star formation that are defined by HIT 
regions in Ha images of other galaxies tend to be concen- 
trations of 0-type stars, and therefore have ages of around 
10 million years. These regions are the classical OB asso- 
ciations. Because of the At — S correlation, they have a 
characteristic size that corresponds to their age. For a mea- 
sured characteristic size of OB associations equal to ^ 80 pc 
(Lucke & Hodge 1970; Efremov, Ivanov, & Nikolov 1987), 
equation ^ confirms that their duration of star formation 
is ~ 14 Myr. The identification of these regions is entirely 
based on the presence of O-stars and bright emission nebu- 
lae, and is therefore only a selection of one particular scale 
out of a continuum of scales for the star-formation process. 
This was implicitly the case in Efremov et al. (1987) and 
Battinelli et al. (1996), where stars in M31 were selected to 
be brightest in U or B to detect the 0-associations. However, 
OB associations are not representative of the star formation 
process in general; they are only one level in a continuous 
hierarchy of self-similar processes that extends from parsec 
to kiloparsec scales (Elmegreen & Efremov 1996, 1998). 

Figure ^ shows a schematic diagram of the sizes and 
durations of star formation in regions that are commonly 
discussed. The At — S relation here is based on the square 
root approximation to the power in the size-linewidth re- 
lation (cf. Eq. I) and not on the slightly shallower power- 
law dependence for clusters written in equation (^. The 
square root dependence gives a more sensible result over 
the whole range of star-forming scales. The figure considers 
small clumps in which individual stars form on 10^ year time 
scales, T-Tauri stellar associations, OB subgroups and asso- 
ciations, and star complexes, in which stars typically form 
on 1, 3, 10 and 30 Myr time scales, respectively. The figure 
also includes spiral arm segments (not density waves) that 
may extend for ~ 10 kpc. The ages of spiral arm segments 
come from the pitch angles (~ 15°) and the rate of shear in 
the galaxy. We discussed previously how such star-forming 



spiral arms, which are clearly distinct from density waves 
by their short lengths and lack of dustlanes, are a natural 
extension of the star formation process to scales larger than 
the disk thickness (Elmegreen & Efremov 1996). 

5. Star formation may be intimately connected with tur- 
bulence because of the similarity between the At — S relation 
and the S/c — S relation for molecular clouds. If true, then 
turbulence would contribute to the rate of star formation in 
a region. The formation of cloud structure by turbulent pro- 
cesses is well recognized (e.g., Langer et al. 1995; Falgarone 
& Phillips 1996). 
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Table 1. Cepheids in LMC Clusters with Integral UBV Photometry 
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